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Low energy spectrum of the S = 1 kagome Heisenberg antiferromagnet (KHAF) is studied 
by means of exact diagonalization and the cluster expansion. The magnitude of the energy 
gap of the magnetic excitation is consistent with the recent experimental observation for m- 
MPYNN-BF4. In contrast to the 5 = 1/2 KHAF, the non-magnetic excitations have finite 
energy gap comparable to the magnetic excitation. As a physical picture of the ground state, 
the hexagon singlet solid state is proposed and verified by variational analysis. 

KEYWORDS: kagome Heisenberg antiferromagnet, Pade approximation, cluster expansion method, numerical 
diagonalization, hexagonal singlet solid state 



§1. Introduction 

The quantum antiferromagnets on the kagome lattice 
has been extensively studied theoretically and experi- 
mentally because of the interest in the interplay of the 
strong quantum fluctuation and the highly frustrated na- 
ture of the lattice structure. So far, most flf-,the attempts 
has been focused on the 5=1/2 caseBaBtf Owing to 
these extensive studies, it is widely accepted that the 
ground state of this model is a spin singlet state and the 
triplet excitation has a finite energy gap. Furthermore, 
there are a number of singlet excitataions below the first 
triplet excitation possiblly down to zero energy in the 

thermodynamic limit .If 
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Fig. 1. The structure of m-MPYNN-BF 4 lattice. 



Recently, Wada and coworkersEmLP have synthesized 
the material m-MPYNN-BF 4 and measured its magnetic 
and thermal properties. This material has the two di- 
mensional structure which is schematically depicted in 
Fig. 1. Magnetically it can be regarded as the Heisen- 
berg model consisting of dimers of two 5=1/2 spins on 
each bond of the triangular lattice. The intradimer ex- 
change coupling Ji is ferromagnetic and depicted by the 

* e-mail: hida@phy.saitama-u.ac.jp 



thick solid lines. These spins are also coupled via inter- 
dimer antiferromagjietic couplings Ji represented by the 
thick dotted lines fltjEP Thus the Hamiltonian is given 
by 



H = Jl £ 



SiS<i 



■h 



E 



(i.i) 



intradimer 



intcrdimer 



where Si is the spin-1/2 operartor. Experimentally, the 
exchange constants are estimated as J\ ~ — 23.26K and 
J 2 — 3.1 IK from susceptibility measurements, so that 
| J\ \»\ J2 |- Therefore, at low temperatures, this 
material can be regarded as the spin-1 kagome Heisen- 
berg antiferromagnet (KHAF) described by the following 
Hamiltonian, 



H = J 



E 

<id> 



Si S 7 



(1.2) 



where Si is the spin-1 operator and the summation is 
taken over all nearest neighbour pairs of sites of the 
kagome lattice. The effective exchange constant between 
5=1 spins is given by J = J 2 /4 ~ 0.78K. Wada et alB 
confirmed that the magnetic excitation of this material 
has a finite energy gap A ~ 0.25K from low temperature 
susceptibility measurement. n 

Theoretically, Asakawa and SuzukiB' have carried out 
the cluster expansion calculation for the ground state 
energy up to the fourth order in A in their pioneering 
work. They calculated, however, the excitation gap only 
up to the first order so that conclusive result for the 
energy gap was not obtained. 

In the present work, we investigate the ground state 
and excitation spectrum of this model (1.2) by the nu- 
merical diagonalization and series expansion method. We 
also propose a possible physical picture of the ground 
state named hexagon singlet solid (HSS) state. To con- 
firm this picture, we perform the variational analysis of 
the ground state energy based on the HSS wave function. 

This paper is organized as follows: In §2, the model 
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Hamiltonian is presented. The numerical method em- 
ployed in this work is explained in §3. The numerical 
results are presented in §4. The results are compared 
with experimental value for m-MPYNN-BF 4 . The HSS 
picture of the ground state is presented in §5. The last 
section is devoted to a summary and discussion. 

§2. Model Hamiltoninan 

Let us rewrite the Hamiltonian of 5 = 1 KHAF as 
follows, 



H — Tin + XHd, 
7~tu = J ^ ^ SiSj., 



Ha 



(2.1) 
(2.2) 

(2.3) 



where ^ A and represent the summation over the 
bonds in the upward pointed triangular clusters (UTC) 
and downward pointed triangular clustrs (DTC), respec- 
tively, which are depicted in Fig. 2. Although we are in- 
terested in the uniform KHAF, we introduce the trimer- 
ization parameter A for the computational purpose. In 
the following, we take the unit J = 1. 
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Fig. 2. Structure of the kagome lattice, 
are denoted by u and d, respectively. 



The UTC's and DTC's 



§3. Numerical Methods 

3.1 Exact Diagonalization 

For N = 12, 15 and 18, the excitation energies of the 
low lying singlet and triplet excited states are calculated 
by the numerical diagonalization method using the Lanz- 
cos algorithm. The clusters shown in Fig. 3 are used with 
periodic boundary condition. 



3.2 Cluster Expansion Method 

As a complementary method, the magnetic excitation 
with 5 = 1 is calculated using the cluster expansion 
method treating 7id as a perturbation. 

Following Asakawa and Suzuki,tf we take A = in the 
unperturbed state. In this case, the ground state is an 
simple assembly of independent singlet states on UTC. 
In the following, this state is called upward pointed tri- 
angular singlet (UTS) state and the corresponding state 
on DTC's is called downward pointed triangular singlet 
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Fig. 3. The clusters used for numerical diagonalization. 



(DTS) state . We apply the series expansion with respect 
to A using the method of Gelfand, Singh and HuseoEP 
and GelfandEr to obtain the series for the triplet exci- 
tation energy e(fe) up to the third order in A where k is 
the wave vector associated with the excitation. 



e(fc;A) 



M 

E 

m=0 



c m (k)\ m withM = 3. (3.1) 



The energy gap evaluated from this series is in itself ill- 
behaved for A>0.3. We therefore employ the Pade ap- 
proximant, 



[K,L] _ J2k=oPk(k)\ 



E£=o«(fc)A J 



(3.2) 



with po(k) = 1. 

In the usual Pade approximation, K + L must be equal 
to M. The original Hamiltonian (2.1), however, has the 
symmetry H(l/\) = H(X)/X. Therefore the energy gap 
has also the same symmetry e(k; 1/A) — e(k;X)/X. We 
therefore i mpos e the same symmetry for the Pade ap- 
proximant (3^2) which ensures pk = PK-k and qi = qL-i 
with K = L + 1 = M. In the following, we call this 
approximation as 'symmetric Pade approximation'. We 
also applied this approximation to the 4-th order series of 
the groundjState energy per site e g obtained by Asakawa 
and Suzuki® which yields the value e g = —1.4468. This 
is close to the values estimated from exact diagonalzation 
of N — 18 clusters, namely e g = —1.4511 (cluster(a)) 
or —1.4393 (cluster(b)). This confirms the reliability of 
symmetric Pade approximation. 

§4. Numerical Results 

The A-dependence of the energy gaps is shown in Fig. 
4. According to the diagonalization results, the ground 
state is always the singlet state with total spin S — 0. 
The A-dependence of the energy gap is common to all 
diagonalization data and series expansion data for small 
A. For small A, the lowest excitation has always 5=1. 
This corresponds to the triplet excitation localized in 
UTC. As A increases, however, the properties of the low- 
est excited states vary depending on the methods used, 
size and form of the clusters. According to the diago- 
nalization data, the singlet excitation with 5 = comes 
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down as A approaches unity. At A = 1, the lowest exci- 
tation has 5 = 1 or 5 = depending on the system size 
and forms of the cluster. They are summarized in table 
I. The first excited state has 5 — 1 for TV = 12 and type 
(a) cluster with N = 18. On the other hand, for N = 15 
and type (b) cluster with N = 18, the first excited state 
is 5 = non-magnetic state. The magnetic gap and non- 
magnetic gap are close to each other for N = 18. As for 
the wave number of the excitations, most of the exact 
diagonalization data shows that the lowest excited state 
has wave number k = (T-point) except for the case of 
type (a) N = 18 cluster for which the magnetic gap is lo- 
cated at the K-point near A = 1. On the other hand, the 
series expansion result shows that the magnetic gap is 
located at the W-point of the Brillouin zone for A > 0.48 
while it is located at the origin for A < 0.48. The dis- 
persion relation estimated from the series expansion is 
shown in Fig. 5. 




Fig. 4. The A dependence of the energy gaps. The small sym- 
bols are the estimation from the symmetric Pade approximant 
of the cluster expansion series. The large symbols are the ex- 
act diagonalization results for N = 12, 15 and 18, respectively. 
The inset is the Brillouin zone of the kagome lattice. Energy 
gaps are normalized by a factor (1 + A)/2 so that the symmetry 
A(l/A) = A(A)/A is evident. 




Fig. 5. The dispersion relation of S = 1 KHAF estimated from 
the symmetric Pade approximant of the cluster expansion series. 
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Table I. Energy gap, spin quantum numbers and wave numbers 
of the low lying excited states of finite size clusters. The values 
with * are the lowest excitations for each cluster. 



N 


A(5 = 0) 


k 


A(5 = 1) 


k 




12 


0.67074 





0.67051* 





15 


0.30050* 





0.48270 





18a 


0.49656 





0.42262* 


K-point 


18b 


0.38134* 





0.42343 






Fig. 6. The system size dependence of the energy gaps for A = 1 
for (a) S = 1 and (b) S = 1/2 KHAF. The small symbol in 
(a) is the estimation from the symmetric Pade approximant of 
the cluster expansion series. The large symbols are the exact 
diagonalization results. In (a), for N = 12 the filled circle is 
almost covered by the open circle. 



To make clear the difference between the excitation 
spectra of 5 = 1/2 and 5 — 1 KHAF, the system size 
dependence of the gaps are shown in Fig. 6(a) and (b) 
for 5 = 1 and 5 = 1/2 cases, respectively. For 5 = 1/2 
KHAF, the non-magnetic gap is much smaller than the 
magnetic gap. Actually, it is expected that the low 
energy non-magnetic excitations are distributed quasi- 
continuously down to zero according to the detailed cal- 
culation with larger clustersE'CF For 5=1 case, how- 
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ever, both gaps have the same order of magnitude within 
the system sizes studied here and we find no indication 
of the appearance of the quasi-continous non-magnetic 
excitation spectrum. 

In any case, the magnitude of the magnetic gap is es- 
timated between 0.42 J (N = 18 cluster) and 0.26J (se- 
ries expansion). The magnetic gap for m-MPYNN-BF^s 
estimated as 0.25K from low temperature susceptibility 
measurement. This value corresponds to 0.32 J which is 
consistent with present estimation. 

§5. Hexagon Singlet Solid Picture 

The presence of finite magnetic gap suggests that the 
ground state consists of an assembly of localized singlet 
states. Apparently, there are three candidates of such 
singlet states. Two of them are UTS and DTS states 
which are the ground states for A = or A — > oo, re- 
spectively. At A = 1, however, the symmetry betwenn 
UTC's and DTC's is recovered. The ground state there- 
fore should contain both states with equal weight. 




Fig. 7. The hexagon singlet solid state. The six S = 1/2 spins 
on a circle form a hexgaon singlet state. The pairs of S = 1/2 
spins connected by the thick solid lines are symmetrized to yield 
5 = 1 spins. 



There is another possibility that the singlet states are 
formed on the hexgaonal clusters. Especially, if we de- 
compose a 5 = 1 spin into two 5 =; 1/2 spins, as in the 
VBS model of Haldane gap state,E3 we can construct a 
trial state covered by local singlets without breaking the 
UTC-DTC symmetry as depicted in Fig. 7. Hereafter 
we call this state as hexagon singlet solid (HSS) state. 

The wave function of the HSS state can be written 
down by expressing the single 5 — 1 spin state as, 

I 1 >= ^1/2,1/2/^2 (5.1) 
I >= $1/2,-1/2 = "0-1/2,1/2 
| -1 >= 0-1/2,-1/2/ 'V2 

where ip a ^ is the basis state of the 5 = 1 spin in terms 
of the two symmetrized S =1/2 states defined as, 

Q ,/3 = ['0a®^/3+'0/3®'0«]/v / 2. (5.2) 

Here ipa is the basis state of the 5=1/2 spin. 

Following the construction of the VBS state by Affleck 
et aljEP we construct the HSS state as, 



$ HSS = (g) H(s auli + * Aj74 ) 

i i 

x w^ ip ' rkp ~ ll p~ <m p (5-3) 
p 

where the summation is assumed for the repeated in- 
dices. The i-th site has two 5 — 1/2 indices on and 
Pi. The index p distingush the hexagons and the sites 
{i v , j v ,kp,lp,mp,np\ belong to the p-th hexagon. The 
wave function yf i p^ip^ k v^ l p^ m v^ ri p is that of the lowest 
singlet state of a single hexagon composed of six 5=1/2 
spins. 

Taking into account the two triangular singlet states 
$ UTS and $ DTS , we take the variational state as 

$ = c u $ UTS + c d $ DTS + c h $ HSS (5.4) 

and minimize the energy expectation values for two types 
of clusters with N — 18 for A = 1. As shown in table 
II, c u , Cd << Ch and the variational ground state energy 
is close to the exact value for both clusters. If we only 
admit two triangular singlet state, the variational energy 
per site is almost equal to — 1 which is far from the exact 
value. On the other hand, if we take only the HSS state 
the variational energy is around —1.3 for both clusters. 
This value is already very close to the optimum value. 
Thus the HHS picture gives a reasonable physical picture 
of the ground state of 5 = 1 KHAF. 

§6. Summary and Discussion 

The spin-1 kagome lattice antiferromagnetic Heisen- 
berg model is studied by means of the cluster expansion 
and exact diagonalization method. The quantum num- 
bers of the lowest excited state depends crucially on the 
method used, size and form of the clusters. This, in 
turn, confirms that the excitation spectrum of KHAF is 
strongly degenerate even for 5 — 1 case. It is suggested 
that the nonmagnetic excitations have finite energy gap 
as well as the magnetic excitations. The HSS picture of 
the ground state is presented and varified by the varia- 
tional argument. 

The material m-MPYNN-BF.4, which motivated the 
present study, can be regarded as the assembly of 5 — 
1/2 spins on the hexagon as shown in Fig. 1. If the 
interhexagon coupling J\ vanishes, the ground state is 
the six spin singlet state on each hexagon. If J\ is 
switched on, the complete singlet state is perturbed. As 
far as J\ is ferromagnetic, however, it cannot kill the 
spin degrees of freedom totally. Actually, we know that 
the 5=1/2 antiferromagnetic-ferromagnetic alternating 
Heisenberg chain has an energy gap even in the limit of 
infinitely strong ferromagnetic bond giving rise, to the 
Haldane phase of the 5=1 Heisenberg chainO From 
this analogy, also, it is reasonable that the magnetic gap 
remains finite in the present material even in the limit 
of J\ — > — oo which is nothing but the 5=1 KHAF. 

On the basis of HSS picture, we can also understand 
the physical origin of the finite non-magnetic gap in the 
following way. There are local singlet excited states in 
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Table II. Various estimations of ground state energy per site e g for A = 1. The values estimated from variational function (5.4), pure 
HSS state, UTS+DTS state, numerical diagonalization and series expansion are tarbulated. For the variational estimation, optimum 
values for the coefficients of HSS, UTS and DTS states are also shown. 



Cluster 


e g (HSS+UTS+DTS) 




Cu(= C d ) 


e g (HSS) 


<E g (UTS+DTS) 


€g (exact) 


e g (series) 




18(a) 


-1.3031 


0.9921 


0.0326 


-1.3028 


-1.0015 


-1.4511 


-1.4468 


18(b) 


-1.3081 


0.9966 


0.0168 


-1.3080 


-1.0010 


-1.4393 





each single S = 1/2 hexagon. These states have the 
exciation energy of the same order of magnitude as the 
magnetic excitaions. If J\ is switched on, these states 
will tend to a gapped non-magnetic excitations. 
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